Quantum Separability of Thermal Spin One Boson Systems. 
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Using the temperature Green's function approach we investigate entanglement between two non- 
interacting spin 1 bosons in thermal equilibrium. We show that, contrary to the fermion case, the 
entanglement is absent in the spin density matrix. Separability is demonstrated using the Peres- 
Horodecki criterion for massless particles such as photons in black body radiation. For massive 
particles, we show that the density matrix can be decomposed with separable states. 
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Recent progress in quantum information theories and experiments 



non-locality and entanglement in many particle systems J, [5|, [6 




has led to interests in studying the 



16j such as Bose 



Einstein condensations 17J, Heisenberg Model[18|, Fermion systems [19|, |20|, and superconductors [2l| and even in the 
vacuum 22| . Entanglement is now treated as a physical quantity like energy and entropy, as well as a resource for 
quantum information processing. Spin 1 particles such as photons and W and Z a gauge vector bosons are essential 
ingredients in the standard model. Furthermore, black body radiation(BBR), the historical birth place of the quantum 



231 ] . the black hole radiation 



2J and the 



physics, still plays an important role in many fields such as quantum optics 
cosmic microwave background radiation [25j. Thus, studying entanglement of thermal spin 1 bosons is important. 

In this letter, we use the temperature Green's function approach to investigate the quantum entanglement of 
two non-interacting (massless and massive ) spin 1 boson particles in thermal equilibrium. Vedral[19j studied the 
entanglement in many body systems at zero temperature using the second quantization formalism. Following his 



works, Oh and Kim 



and the Kondo model 



20] studied the entanglement of two electron spins in a free electron gas, superconductivity 



21] 



261 ] at finite temperature using thermal Green's function methods. In his work, Vedral showed 



that there is no reason why the polarizations of a pair of separated photons should be correlated at zero temperature. 
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At finite temperature, however, the situation becomes more complicated. In this case entanglement may occur because, 
contrary to the intuition that thermal noise destroys entanglement, it has been shown that even if two particles do not 
interact directly, they can become entangled by interacting with a common heat bath 



27 



28 1 . To maintain thermal 



equilibrium, the thermal bosons (like photons in the BBR) should interact with a common thermal bath, even when 

they are not directly interacting with each other. Therefore we need to calculate the entanglement of thermal spin 1 

bosons explicitly to check for the absence of entanglement in the systems. Recently, it was also shown that two qubits 

□ 

interacting with the BBR can be ent angled [29J. In earlier works entanglement in many body systems was usually 
tested indirectly by investigating the entanglement of two'probe qubits' interacting with the system. In this letter, 
however, we are interested in the entanglement of the particles themselves without any probe qubit. This approach 
could reveal the physical nature of the system more clearly. 

We begin by briefly reviewing the Green's function approach 30j • To calculate the entanglement we need to know 



the density matrix of the system with Hamiltonian H and temperature 1//3. The finite temperature two-particle 
density matrix is defined with the field operator ip(xi) for the i-th particle; 

pW{x u x 2 ;x' x x' 2 ) = ^Va)^Vi)fe)to)> (1) 

= --y{X 1 T 1 ,X2T2;X 1 T{ r ,X 2 T 2 r ) 

where (O) — Tr{paO} with Z = Tr{e~' 3ff } and pa = e~^ H /Z, and r/~ (i = 1,2) denotes a time infinitesimally later 
than t,;. Using the Wick's theorem, the two-particle temperature Green's function for bosons can be reduced to the 
product of one-particle Green's functions; 

^(W, 2')= Tr {p G T T [^(l)^(2)^(2')^(l')]} (2) 
« 0(1, l')S(2, 20 + 0(1, 200(2, 1') 

where the number % (i = 1,2) denotes the space-time coordinates (0^,7$) of particle i, and 0(1; 1') = 
Tr{/5c T T [4>k( 1)V , ^-(1')] } i s the one-particle temperature Green's function. The field operator is redefined as 
$k{x,t) = e Kr l h ip(x) e~ KT l h with K = H — pN, where p is the chemical potential and N is the number oper- 
ator. The second equality of Eq. ([2]) denotes the Hartree-Fock approximation which is exact for non-interacting 
systems such as the one considered in this paper. Then, the non-interacting one particle Green's function 0°(1; 1') is 

P {1) (x; x') = -G°(xr; x't+) = <W 9 (r - O, (3) 



where a denotes the spin index and g(r — r') is the one-particle space density matrix in a volume V; 

5( r ) = ^E eIkr ^- (4) 



V 

k 



Here = {exp[/3(ek — (j)] — 1} 1 is the mean occupation number in the state with momentum k and energy 
£k = h 2 k 2 /2m for massive non-relativistic bosons or £k = Kk/c for photons. With Eqs. {2} and ((3|), one has the 



explicit form for the two-particle space-spin density matrix 



31 



32] 



p {2) (x 1 ,x 2 ;x' 1 ,x' 2 ) = hg{ri - r'Jgfa - r'^S^S^ + g(ri - r' 2 )g(r 2 - r'^S^S^} . (5) 

where <Ji denotes the spin index for the z-th particle. To the best of our knowledge, there is still no consensus on 
how to deal with entanglement between continuous variables such as coordinates and discrete variables such as spin. 
Hence we set ri = and v 2 = r' 2 to consider only discrete (spin) degrees of freedom, which leads to a simpler form 
for the density matrix. For isotropic cases, the two-spin density matrix, depending on the relative distance between 
two particles r = |ri — r^l, is 

where a is the number of spin degrees of the freedom (a = 2 for massless spin 1 bosons and 3 for massive ones) . n = 
N/V is the particle density for particle number N and f(r) is an exchange term representing the indistinguishability 
of bosons: 

k 

A bipartite state p is called separable if it can be written in the form 

n 

P = y £,Pipf®pf, (8) 

i=l 



where pf and pf are states of subsystem A and B, respectively. We use the Peres-Horodecki separability criterion 



33, 



34j, which is the positive partial transposc(PPT) criterion. A state is PPT if p TB > 0, where the partial transposition 



of p is 

Pjmjn = (h m \p TB \j, n) = p in ,j m (9) 

in some basis. Let us first consider massless particles such as photons from BBR with two spin degrees of freedom 
denoted by two level states(|0), |1)). The two-spin density matrix corresponds to the two qubits density matrix in this 



4 

case. By dividing the bracket part of Eq. © by 4 + 2/ 2 , we obtain the normalized two-spin density matrix p\2 32 1 
for a given relative distance r between two photons in {|00), |01), |10), |11)} polarization basis 



1 



Pl2 = 



4 + 2/ 2 



1 + / 2 
1 P 
p 1 



(10) 



o i + p 

where Tr,j l0 . 2 {pi2} = 1 and we have dropped r in f(r) for simplicity. This matrix has the same form as for the 



20(. One can easily show that pi2 is PPT and hence 



Fermion case except that the off-diagonal terms have plus signs 
separable 35 1 (The lowest eigenvalue of pjf is 1/(4 + 2p) > 0). Hence we can conclude that there is no entanglement 
in the two-spin density matrix of the non-interacting massless thermal spin 1 boson system. 

What can the absence of entanglement in the BBR be used for from a practical viewpoint? The absence of 
quantum correlation can help us to understand the nature of certain light sources, for example, astronomical objects. 



By performing an Aspect-type Bell test experiment 

n 



361 ] on polarization states of two photons from a light source and 
checking for violation of the Bell inequality [37|, one could determine whether the source emits entangled photons. If 
this test reveals pairs of entangled photons from the source, one can say that the source is, at least, not a black-body 
radiator. Given that information from distant astronomical objects is mainly obtained by observing electromagnetic 
waves, this quantum test would provide us additional useful information about the objects. 

We now move on to the case of massive spin 1 particles, which have 3 spin states (a = 3). In this case Eq.© reads, 
in {|00),|01),|02),|10),---,|22)} basis, 



Pl2 



9 + 3/ 2 
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(11) 



, which describes a two-qutrit mixed state. Now p^f has the smallest eigenvalue 1/(9 + 3/ 2 ) > 0,and hence pi2 is 



PPT. Since PPT is a necessary, but not a sufficient, condtion[38| for separability in a 3x3 system, PPT alone does 
not automatically guarantee separability. We therefore need to show explicitly that there is no bound entanglement 
in pi 2 . Indeed, it can be shown that p i2 can be decomposed with separable states (i. e., in the form of Eq.®) as 
follows; 



where 



1 



Pl2 = 



9 + 3/ ; 



:(9/ 2 Po + 3o-o + 6(1-/ 2 )cti), 



(12) 



a = -(|00)(00| + |11)(11| + |22}(22|) and 

a x = -(|01)(01| + |02)(02| + |10)(10| + |12)(12| + |20)(20| + |21)(21|) 
are trivially separable. Exploiting roots of unity, po can be written as an integral over product states; 



(13) 



2 * 

Z7T 



(14) 



where \i/>(0)) = ^(|0) 



ID 



!2». 



n 



The massive spin 1 boson class contains vector mesons and massive gauge bosons such as W and Z [39[, which 
intermediate forces between subatomic particles. They are believed to have been in thermal states in the early 
universe. Another interesting example is spinor Bose Einstein condensation with alkalis with hyperfine spin 1 such 
as 23 N a and 39 K 



4fJ. Since these particles usually have self-interactions, it might be interesting to study how 
self-interactions change our results (Our results correspond to tree- level approximations in these cases). 

It is well known that, in algebraic quantum field theory, due to the Reeh-Schlieder theorem[41] the vacuum for 



quantum fields could have quantum nonlocality 



42, 



43j entanglements of Bose (and Fermi) 



43j . Furthermore, in Ref . 

quasi-free (gaussian) states are extensively studied. Our results do not contradict these results, because in our model 
we consider only the discrete (spin) degree of freedom of two separated particles in a thermal system. In our approach 
separability of the two-spin density matrix is shown explicitly using the familiar Green's function method, especially 
for the 3x3 thermal vector boson system. 

Several points can be made in relation to the material presented here. First, using the Green's function method 
we showed that there is no entanglement in the two-spin density matrix of a non-interacting thermal spin 1 bosons, 
regardless of their masses. However our results do not rule out the possibility of more sophisticated entanglements 
in these systems, such as those between continuous variables and discrete variables 
interacting with BBR become entangled 29], even though the spin states of two photons in BBR are separable? 



44|. How do two probe qubits 



One possible explanation would be that BBR has a more sophisticated entanglement, as described above. Another 
explanation would be that the coupling between two probe qubits and the common environment (like BBR) induces 



an indirect interaction between the pair of qubits 



29(. Second, the absence of 



45j . which produces the entanglement 
entanglement of two photon spins can be useful for characterization of light sources by checking the validity of the 
Bell inequality experimentally. Finally, it will be interesting to investigate how self-interactions or higher spins of 
bosons changes the results. 

P.S. After submission of our paper we found that there appears a paper about spatial entanglement of free thermal 



bosonic fields (quant-ph/0607069) 
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